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Abstract 

We investigate the generation of large scale magnetic fields in the universe 
from quantum fluctuations produced in the inflationary stage. By coupling these 
quantum fluctuations to the dilaton field and Ricci scalar, we show that the mag- 
netic fields with the strength observed today can be produced. We consider two 
situations: First, the evolution of dilaton ends at the onset of the reheating stage. 
Second, the dilaton continues its evolution after reheating and then decays. In 
both cases, we come back to the usual Maxwell equations after inflation and then 
calculate present magnetic fields. 
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1 Introduction 



It is well understood that magnetic fields are present on various scales in the universe. 
These scales vary from planet size to clusters of galaxies of Mpc order[l]. The problem 
of obtaining a reasonable mechanism for generation of such a large scale magnetic field 
is an important question in modern cosmology because of its direct influence on the evo- 
lution of universe and astrophysical events. Magnetic fields have an important role in 
the dynamics of galaxies by confining the cosmic rays or transferring angular momentum 
away from protostellar clouds so that they may collapse to become stars (without the 
loss of angular momentum, protostellar clouds would collapse to a low-density centrifu- 
gally supported, unstarlikc state). They also play an important role in the dynamics of 
pulsars, white dwarfs, and even black holes. The strength of these magnetic fields varies 
from /iG in galaxies and cluster of galaxies to a few G in planets and up to 10^^ G in 
neutron stars. 

Some mechanisms known as galactic dynamo are established to amplify the scale and 
strength of magnetic fields by transforming the kinetic energy of turbulent motion of 
interstellar medium into magnetic energy. The gravitational collapse of matter to form 
galaxies and cluster of galaxies is definitely an amplification mechanism that can affect 
the strength of these fields because of the magnetic flux conservation. These mecha- 
nisms, however, are only amplification mechanisms and require a seed magnetic field to 
amplify to strengths presently observed . Theories established for generation of these 
seed fields are classified into two groups: astrophysical processes and cosmological pro- 
cesses in the early universe. 

If the scale of magnetic fields of galaxies and cluster of galaxies are in the order of Kpc 
and Mpc, it means that we should investigate their origin in the early universe rather 
than in astrophysical processes. Then, these seed magnetic fields are trapped in highly 
conductive plasma collapsed to form structures like galaxies and their clusters during an 
adiabatic compression and, finally, subjected to some amplification mechanism like the 
galactic dynamo. 

Many different mechanisms have been proposed for generation of these seed magnetic 
fields in the early universe [4] that may be classified as follows : 

1. Breaking of conformal invariance of electromagnetic interaction at inflationary 
stage. This can be realized either through new non-minimal (and possibly non gauge 
invariant) coupling of electromagnetic field to curvature [5], or in dilaton electrodynam- 
ics [6], or by the well known conformal anomaly in the trace of the stress tensor induced 
by quantum corrections to Maxwell electrodynamics [7]. 

2. First order phase transitions in the early universe producing bubbles of new phase 
inside the old one [8]. A different mechanism but also related to phase transitions is 
connected with topological defects, in particular, cosmic strings [9]. 
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3. Creation of stochastic inhomogeneities in cosmological charge asymmetry, either 
electric [10], or e.g. leptonic [11], at large scales which produce turbulent electric cur- 
rents and, in turn, magnetic fields. 

It seems that inflation is the most natural way for overcoming the large correlation 
scale [5, 12]. Inflation produces effects on scales much larger than Hubble horizon in a 
natural way. Then, if electromagnetic quantum fluctuations had been produced in that 
epoch, they could have been present as large scale magnetic fields today. The idea is 
based on the assumption that a quantum mechanical mode (in scales much smaller than 
Hubble Horizon) is excited which freezes when passing through the horizon. The prob- 
lem arises from the fact that conformally invariant theory can not produce nonmassive 
particles in a conformally fiat gravitational background. This is what happens to pho- 
tons in a FRW background and, therefore, electromagnetic fields can not be produced. 
And if the origin of magnetic fields of galaxies and cluster of galaxies were quantum 
fluctuations produced in the inflationary stage, then conformal invariance of Maxwell 
theory should have been broken in that epoch. This happens in several ways [4]. 

We break the conformal invariance of Maxwell theory by coupling gravity (Ricci 
scalar) and a scalar field (dilaton) to it. A non-minimal coupling of electromagnetic 
fields to gravity was introduced first by Turner and Widrow [5]. Also, the coupling of a 
scalar field to electromagnetic fields has been studied by different people [5, 6, 13, 17]. 
We mix these two situations and consider the more realistic condition that both gravity 
and a scalar field, that is dilaton field, are coupled to electromagnetic field during the 
infiationary era. In [5], the problem is solved qualitatively and different values for mag- 
netic fields are derived by changing the parameters. On the other hand, the problem is 
treated parametrically and more quantitatively in [13]. By integrating these solutions, 
we get a more generalized situation in which the coupling parameters are fixed by using 
some special values for the magnetic field. 

We first introduce our theory, then derive equations of motion and solve them to get 
an expression for electromagnetic vector potential, A^. Then we derive the evolution of 
electric and magnetic fields before and after inflation by using joining conditions. We 
assume two situations here : 

1. Evolution of dilaton field ends by finishing infiation when dilaton freezes. 

2. A more realistic situation in which dilaton continues its evolution after reheating 
and then decays into radiation. 

Finally, we calculate today's strength of magnetic fields on different scales. 
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2 Action and the equations of motion 



For investigating the evolution of magnetic fields produced by a Maxwell theory whose 
conformal invariance is broken, we first introduce the lagrangian. We will then derive 
the equations of motion [5, 6, 13]. 

We introduce inflaton and dilaton scalar fields lagrangian density as 

jC,u^-lg''''d^^d,^-vm , (1) 

An/ = -^^'^'^9^9,0 (2) 

V[^] = Vexp{-XK^) . (3) 
C/[0] and V[^] are infiaton and dilaton potentials, A is constant and k, — The 

■'"pi, 

forms of f/[0] and are determined from higher dimension theories reduced to four 
dimensions [15, 16, 17]. 

The lagrangian density of electromagnetic field coupled to the dilaton and gravity is 

Cem = fmi—F.^F^'' + iRA^) , (4) 

i^=6( '"+;J + ^ ) . (6) 

/($) and ^RA^ are dilaton and gravitational couplings to electromagnetic field, k is 
curvature constant and R is Ricci scalar [6]. a is scale factor and ^ is a dimensionless 
parameter which will be determined. 

Therefore, the action is 

5 = y d^x^/^ [Cinf + Cdu + ^Cbm] ■ (7) 

where g is the metric tensor. We consider the metric of space-time as fiat FRW metric 
(k=0) 

rfs^ = gf.udx'^dx'' = dt^ - a^it) dx^ . (8) 
Our work is based on four assumptions: 

1. In the infiationary stage with slow roll condition, the energy density of infiaton 
field is much bigger than that of the dilaton field, P0 ^ p$. 

2. The universe becomes immediately hot after the infiationary stage, t > Ir. 
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3. The conductivity of the universe is ignorable in the inflationary stage because 
density of charged particles is very small in that epoch. After reheating, a lot of charged 
particles are produced so that conductivity immediately jumps to a large value after 
inflation, Uc S> H. 

4. We consider two different situations for dilaton. First, the evolution of dilaton 
fleld ends by finishing infiation and the dilaton freezes; therefore, the coupling will be 
removed (/ = 1) [6]. Second, the dilaton continues its evolution after reheating until it 
reaches its minimum potential and then decays into radiation when again / = 1 [13]. 



2.1 Equations of motion 

From action (7), the inflaton, dilaton, and electromagnetic potential equations of motion 
are given by 

^ d,{V^99'^dA) + ^^-^. (9) 



'=^iV^99-3.^) + ^ = ^ + ^RA^) , (10) 

^ d,{^f{^)F^'' -2iRA'') = Q. (11) 



Since every inhomogeneity of space will be diluted in the inflationary stage, we can 
ignore the spatial dependence of the inflaton and the dilaton scalar flelds. The right 
hand side of Eq.(lO) is very small and can be considered as a perturbation in the dilaton 
theory. 

The inflaton and dilaton equations of motion are derived from Eqs.(9) and (10) 

^ + 3ff^+!^ = 0, (12) 
$ + 3//$+ ^^^^ = . (13) 
H is the Hubble constant and is derived from Friedman equation 

where 

P4> = + U[ct>] , (15) 



2 

1 . 



= + vm . (16) 
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are the inflaton and dilaton energy densities. Dot implies time derivatives and p — Ptp+p^ 
is total density. According to our assumptions, p^ ^ p$ and we can ignore the dilaton 
energy density in the inflationary stage whereby we will have: 

^ = . (17) 

We have used slow roll condition in the above equation and Hinf is the value of 
Hubble constant in the inflationary stage. Since we assume that the term ^RA^ is small 
compared to F^i,F^^ , we can use the Coulomb gauge, /lo(t, a;) = and djA'it, x) = 0, as 
used in the standard Maxwell theory. We derive equation of motion for electromagnetic 
potential A (in comoving coordinate) from Eq.(ll) 

Mt, x)+(^H+ I j Mt, x) - {^^djdj - 2^R^ Ai(t, f) = . (18) 
2.2 Electromagnetic potential A 

To solve the equation of motion for A, we first quantize the theory. The corresponding 
momentum from jCem is 

T^v — = 1 

dA, 

TTo^O, TT, = /($)a(t)A(t,f) . (19) 



Commutation relation for Aj and tt,- is 



[A^it^xUAt^y)] = ^ J ^se^'-^'-y^ (5., - Mlj . (20) 

where k is comoving wave number. Using these relations, we can write the quantum 
form of A as [17] 



A{t,x) = / ^!^[6(^)A,(i,^)e^*;- + 6t(^)^,*(^^^)e-^-] . (21) 



in which h and 6^ are annihilation and creation operators with the following commutation 
relations 

[h{k),h\k')\ ^5\k-k') , 

[bik);bik')]^[bHk),bHk')]-0. (22) 

For simplicity, we choose xi along the direction of k. Thus, from now on we will work 
only with two transverse components (1=2,3). Equation of motion for electromagnetic 
potential Fourier modes will be: 



Ai(t, k) + ( Hi^f + ^] Ai(t, k)+('^- 2^R ) Ai(t, k)^0. (23) 
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Ricci scalar in the inflationary stage is R — 12H'^ [12]. 

We get normalization condition for Ai{t, k) from Eq.(20) as 

A,{t, k)A*{t, k) - fc) = [si, - ^) . (24) 

For more simplicity in solving the equations, we use the following approximation for the 
evolution of /($) 

fm = nm] = mnm ^ ha'^-' ■ (25) 

where /o is a constant and all time variations of /($) have been put in a^"~^. It can 
be shown that time dependence of a can be ignored if we have slow roll condition for 
inflaton and dilaton fields [13]. The only problem that remains is determination of an 
acceptable range for variation of a that can be achieved from consistency conditions 
(section 4). 



Prom Eq.(25) 



Hinf + ^ = (2a - l)Hir,f . (26) 



By introducing a new variable dr] = ^, and relying on the fact that r] = — in the 
inflationary stage ( in De'sitter spacetime ), we can rewrite eq.(23) as 



One of the forms of Bessel Equation is 

where i^, u' and 5 arc positive and the solution is 

u = z''{DiHI]\Sz) + D2HIJ\Sz)) . (29) 

and Hi)^ are Hankel functions of first and second order u', respectively. Comparing 
Eqs. (27) and (28), we will have 

r] — z V = a 5 — ±k , (30) 



iy'= ± - 24e = ± V"' - 24^ • (31) 
The solution obtained for Eq.(27) is 

Ai{k,rj) = D,{a){-Hir,frjYHi}\-krj) + D2{a){-Hi,,frjyHiJ\-krj) . (32) 
Di and D2 coefficients are determined from the renormalization relation (24) 
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For simplicity, we choose 



Di(a)= /-^/(a) a^-V^e^C'^'+iW^ D,{a) = . (34) 

Since we are working with large scale magnetic fields, we expand Eq.(32) whithin 
large wavelength hmit ior u > and u < 



I ^ e.(.'-iW4 j,/>o, (35) 



aH, 



Aj{k,a) = 2-^\—-^—^V{-y')a-'l^ 

f^— Ve^^^-'^'W^ iy'<0. (36) 

V in/ / 

The electric and magnetic fields can be easily derived from the above expressions 
using the relation between A and E and B in comoving coordinate system. 

3 Electric and magnetic fields evolution 

In Coloumb gauge, electric and magnetic fields are derived from vector potential A as 
follows 

E^—^A, B^VxA. (37) 

In comoving coordinate, electric and magnetic fields are [17] 

Ei^{t,x)^Mt,x) , (38) 



5,^(t,f) = -e,,kdjAk{t,x) . (39) 
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a 



where "c" implies comoving and dj is differentiation in comoving coordinate. 

Using relations (35), (38), and (39), we can derive Fourier components of electric and 
magnetic fields in the infiationary stage 



Ei'^ik, a) = ,/^^^ (-)a-^/^[(- - l/ml'' ' — 



AHinff{a) a " \aHinf ) 



Bi''{k,a)^-i{-l)\ 



TT 



AHinffia) 



X ( -] a-V^i/W ( e^(-'+i)-/4 . (41) 
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We assumed that the dilaton freezes after inflation. In this epoch, the conductiv- 
ity of the universe increases immediately so that a ^ H [6]. And the evolution of 
electromagnetic vector potential follows from the equation 

f) + (- + ae) i.(t, x) - ^djd,A,{t, f) = . (42) 

Ratra has shown that we have the following joining conditions for electric and mag- 
netic fields at transition from the inflationary stage to the radiation dominated epoch ( 
t = tR) [6, 18], 

E'^iiRD) {tR, x) = exp {-actR) E^h^inf) {tR , x) , (43) 

B?^RD) {tR: X) = S,^,jvF) (tR, x) . (44) 

According to these relations, in a universe with large conductivity, the electric field 
accelerates charged particles but it reduces exponentially. Prom Alf'ven theorem, the 
magnetic fiux is conserved in a conductive universe and, therefore, the magnetic field 
evolution is by in the physical coordinate system [6, 17, 18]. 

In the physical coordinate, the electric and magnetic fields are [17] 



Ei''''{t,x)=a-'Ef'{t,x) , (45) 
Bi^\t, x) = a-^Bf'it, x) . (46) 
Using relations (39), (41) and (46) we have 
\Bi''\t,k)\'' = a-^\Bi^{k,a)\'' 

^ {iH^Iffia)) [1] ^^'^ (^) ^^'^ (^) • ^^^^ 

We are considering the magnetic fields on large scales. Thus, we expand relation (47) 
for small k to note that the magnetic field evolves with after infiation 

')2v'-2 

\B[\t,k)\'^- r\u')f-\aR) 



TT 

( \ ( k V^" fi 



anHinf J \aRHinf J \a 



I/' > , (48) 



\B^%k)\' = ^—^T'{-u')f-\an) 

TT 



/ A;^ W k V 

ydRHinf J \aRHinf 



I Q)', u'<0. (49) 
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The magnetic energy density in the Fourier space is 

PB{t,k)^l\B/%k)\' f{a). (50) 

Since two transverse components are equal, Bt = \/2Bi. Multiplying (50) by phase 
space density , the large scale magnetic field energy density in the coordinate space 
is 

ps{L,t) = -^\Bj''\t,k)\'f{a). (51) 
where L = ^ is correlation length. 

By integrating Eqs.(48), (49) and (51), we obtain the magnetic field energy density 
on the scale L — ^ 

,,^,,^^^l_rW>./{^) {^) . (52) 

Accoring to the above relation, the spectrum of the magnetic fields is invariant for 
2\u'\ = 5. 



4 Consistency condition 

The energy density of electric and magnetic fields should be smaller than the energy 
density of dilaton field in the infiationary stage so that only V[^] determines the evo- 
lution of dilaton [17]. We show the ratio of electromagnetic energy density to dilaton 
energy density by Q{L,t): 

e{L,t) = —[pB{L,t)+pE{L,t)], (53) 
PE{L,t) = —\Ej^'{t,k)\'f{a). (54) 

If © ^ 1, the consistency condition is satisfied. Using Eqs.(17), (40), (41), (53) and 
(54), we have 



9w\MpiJ \aHinf J " \aHinf J " \aHinf ^ 



k \ ^^n^ I k \\ ( 1\ ^tA\ ( k 



+ I I 1/2 I Hi;'_, {^\-Hi}U^]] + (u-^] 



(55) 
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As a measure of ©, one may use the upper limit of Hinf calculated from CMB 
anisotropy observations [19, 20, 21] 

< 2 X 10"^ =^ Hir^f = 2.4 X W^GeV . (56) 

Mpi 

and noting that -rr — decreases with evolution of universe, from Eq.(55), we may con- 

elude that when — = 1, © < 10~^^uj~^ ( u is defined as a; = ^ and according 
to our assumption that the dilaton energy density is much smaller than that of inflaton 
in the inflationary era, a; -C 1 ), the consistency condition exists. By investigating these 
relations within large wavelength limits, we have 

97r^w\MpiJ ^ V ; / V V ; 

_22.'+lp^^/ _ ^^^^^f ^^^^ 2^2u'+l)^2^^, ^ ^ 

(i/ - l/2f2'"''T^{iy') + (i/ - l/2)(-22("'+i)r(i/' + l)T{iy')x 

+22'^'r(i/')r(i^' - l)x-')]] iy'>l, (57) 

e ^ -^(EM.Vu-5+^-'2^-'r^(u') + i/4(22(i--')r2(i _ u')x^i^-^^'^ 

-2^r{iy' + i)r(i - 1/')^'^'""' + 2(2-'+i)r2(i/' + 1)^^) + i/2f2'^'T\v') 
+{v - i/2)(22r(i/')r(i - i/')^^'"'"' - 2''^'+^r(i/')r(i/' + i)x)] 

< I/' < 1 , (58) 

^ 1 /':[^')'r3,-(5+2.')2-2-'r2(_^,') + IU(2''^^--')V\1 - u')x^ 
97r^w\MpiJ ^ V ; / V v ; 

-2^r{u' + i)r(i - iy')x^ + 2(2-'+i)r2(j/' + i)x'(i+'^')) + (i/ - i/2)22-2'^'r2(-i/') 

+ (z/ - l/2)(-22r(-z/')r(l + z/')a;i+'"' + 2^^^-'''^T{-iy')T{l - iy')x)] 

- 1< z/' < , (59) 

e ^f:^yb-(5+2.')2-2-'r2(_i.') + i/4(22(--'-i)r2(-i/' - i)x-^ 

_2-2-'+ir(-z/' - i)r(-z/' + 1) + 2(-2'^'+i)r2(-z/' + i),t2) + (z/ - 1 72)22-2^' r2(-z/') 
+{u - i/2)(-22(-'^'+i)r(-i/' + i)r(-i/')x + 2-2^'r(-i/')r(-i/' - 1)^-^)] 

I/' < -1 , (60) 

which wc defined as x = "'^^."■^ . Since ^ decreases with the expansion of the universe, 
a negative sign must be selected for the power of x so that the consistency condition 
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(© <S 1) is maintained. The intervals < v' < 1 and — 1 < i/' < give us trivial results. 

Therefore, v' > 1 and v' < —1 conditions will determine the range of u' by considering 
the consistency condition. From u' > 1 and z/' < — 1, we conclude 2z/' < 3 and 2z/' > —3, 
respectively. By integrating the two situations, the consistency condition implies the 
range —3 < 2u' < 3 for u'. Note that the first terms in Eqs.(57)-(60) is related to the 
magnetic field and the other terms are related to the electric field. 

5 Present magnetic field 

As explained before, the seed magnetic fields produced in the early universe have been 
affected by amplification mechanisms to get the strengths and scales observed today. 
These amplification mechanisms are galactic dynamo and gravitational collapse. 

The galactic dynamo mechanism is based on the conversion of the kinetic energy asso- 
ciated with the differential rotation of galaxies into the magnetic field energy [16, 23, 24]. 
In the ideal situation, this mechanism could have amplified the magnetic field strength 
with a factor of e'^° or 10^^ during the 30 revolutions of protogalaxies since their forma- 
tion up until now [16]. In the real situation, the dynamo mechanism could have amplified 
the magnetic field with a factor of 10^ although it has not been properly obtained in the 
cluster of galaxies. 

The gravitational collapse is another means for amplifying the magnetic fields. For 
the galaxies, right before its formation, a patch of matter of roughly 1 Mpc scale col- 
lapses by gravitational instability. Right before the collapse, the mean energy density of 
the patch stored in matter is of the order of the critical density of the Universe. Right 
after collapse the mean matter density of the protogalaxy is, approximately, six orders of 
magnitude larger than the critical density. Since the physical size of the patch decreases 
from 1 Mpc to 30 Kpc, the magnetic field increases because of the flux conservation, of 
a factor (pa/Pb)^''^ ~ 10^ where Pa ^-nd pb are, respectively, the energy densities right 
after and right before the gravitational collapse [16]. For the cluster of galaxies, their 
density is greater than the critical density of the universe by a factor of 10^ and we have, 

(Pa/Pb)'/=' - 10'. 

The presently observed strengths are 10~^ G in galaxies and 10~^ in cluster of galax- 
ies. By both mechanisms available, one has to have the magnetic fields with strengths of 
10^^"^ G (in ideal galactic dynamo) and lO"^'' (in the more realistic galactic dynamo) for 
galaxies on 1 Mpc scale. When just the gravitational collapse acts on them, we should 
have the seed magnetic fields with strengths of 10~^° G in galaxies and 10~^ G in cluster 
of galaxies (on 10 Mpc scale). 

By noting that in Eq.(52), none of the amplification mechanisms has been considered 
for the current magnetic fields, we consider the values mentioned in the above paragraph 
as the present strengths of the magnetic fields of galaxy and cluster of galaxies and 
emphasize that the observed fields can be achieved by the effect of the amplification 
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mechanisms. 



We get the value of magnetic field at the present time as 



PB{L,to) « \B{L,to)\' « {^Hi^fMpilf'^'^'' ( T,J^] \ (61) 



PI 



where we have used the following relations [12] 



H^lMxgl/' (62) 



PL 



P0 = 30^*^« (9* ~ 200) , (63) 
TV = 45 + ln(^ + m [30/(^V)]^^"p/^^ (64) 



^0 r 9* \1/3^£r_ 



an '3.9r T^o 2.35 x 10-^'^[GeV] 



(65) 



"0" implies quantities at present. T/,> and T^^ are temperatures of CMB at the reheating 
epoch and at the present time, respectively ( T^^ a; 2.73-ftr ). is e-folding between 
time of first crossing, ti, and reheating. We can rewrite eq.(61) as an equality 



i X 4.8 X 10-^9 B{L,tor = ^ ^'^I^'I^sIt)''^' ^"^^ ^"^"^^'^^ 

[77 ^ f_30_N-l/6 l2i/'|-5 

^ Vli^ ) ^ ( ( 10575^^^^^^ '''' ■ '''' 

To find a range for u' variations, the graph of log(i/j„/) with respect to is plotted 
(fig. ??) . By taking the logarithm of both sides, one has 

Log{H,^f) = (^^^) (-(16.041(2 u' - 5) 

+2 Lo5(r(z/') + 0.6 u') + 51.98 - 2Log\B\) . (67) 

As mentioned before, the upper limit for Hinf is 2.4 x lO^^^Gel^. We can estimate a 
lower limit for Hi^j by noting that duration of the inflationary stage had been about 
10^'^'' [22] and that the least value of e-folding number should be about 60 

^ = e^ = e^*^ ^ Hir,f^{Aty'N^AxlO''GeV. (68) 

where "i" means value of quantities at the beginning of inflation. Now, we use Eq.(67) 
to depict the graph of log{Hinf) versus for the magnetic flelds with strengths 10^^'^, 
10"^'', and 10~^° G on 1 Mpc scale (galaxies) and 10~^ G on 10 Mpc scale (cluster of 
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galaxies). Prom Fig. (1), we find that for the magnetic fields mentioned, with 10^^ < 
Hinf < 10^^ the ranges of u' variations are 1.96 < u' < 2.9, 2.21 < u' < 2.37, 2.5 < 
z/' < 2.68 and 2.54 < u' < 2.75, respectively. All of these ranges violate the consistency 
condition, |z/'| < |. From Eqs.(57)-(60), it can be seen that the consistency condition 

comes from H^\^i-^ — ) (which is part of the electric field). For solving this problem, 

we have to cancel out the term Imposing the condition u — ^ = u' on eq.(40) and 

considering the following relations between derivations of Hankel functions 



^i'li(^)+^^ai(^)=2^^«(x) 

X 

d 

dx 



Hi]U^)-Hl]Ux) = 2-Hi]\x). (69) 



the electric field can be written as 



Now rewriting eqs.(57)-(60) using the above relation, we arrive at 

e» ~ 7r^( ^ {x-'+^'''{2^'''r\u') + 22(-'-i)r2(i.' - l)x-2) u' > 1 (71) 



e 



e ^ ^ ^-(5+2.')(2-2-'r2(_i.') + 2-2('^'+i)r2(-i/' - 1)^-^) u' < -^7.3) 



1 ( 






kMpi 


1 ( 


^ Hinf 




yMpi 


1 ( 


''Hinf 




yMpi 



The new consistency condition is ^ < z/' < |. The problem of satisfying the 
consistency condition still remains for u' > ^ and, therefore, the acceptable range for 
u' variation is z/' < |. The new consistency condition is derived by the redefinition of 
i/ — I = i/' while we have a = from Eq.(31) which results in a new consistency condition 
as q; < 3. Using u' — a — ^ in Eqs.(59) and (67) one has 

-1 o2a— 4 —1/3 

i X 4.8 X 10-=^^ BiL,tor = ^ - l/2)(^) x 1.66 T,,iGeV) 

Log{Hinf) = {a- l)-'(-(16.041(2 a - 6) + 2 Log{r{a - 1/2)) 

+ 0.6( a - 1/2)) + 51.98 - 2Log\B\) . (75) 

We conclude from Eq.(75) that the larger Hinf and a, the larger will the magnetic 
field be. For a = 1, magnetic field is 2.1 x 10^^^ G in 1 Mpc scale which is independent 
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of Hinf. For q; = 3, magnetic field spectrum is invariant, and the maximum value of the 
field is 1.8 X 10"^^ G for Hir,f = 2.4 x lO^^GeV. By using Eq.(75) to depict the graph 
of \og{Hinf) versus |z/'|, for the magnetic fields with strengths 10^^^, lO^^'' and 10^^° G 
on 1 Mpc scale (galaxies) and 10~^ G on 10 Mpc scale (cluster of galaxies), Fig. (2), 
it is clear that for 10~^^ G and 10~^^ G fields, the consistency condition holds in the 
appropriate range for Hi^f while it does not for 10~^° G and 10~^ G fields. This should 
not concern us since we got 10~^° G and 10~^ G fields by assuming that the galactic 
dynamo does not act on galactic scales which is not, indeed, realistic. 

As mentioned before, the most real situation for galaxies is a 10~^''G magnetic field 
on 1 Mpc scale, which transforms to a 10~^G field observed today via gravitational 
collapse and dynamo action. In this way, a range for a parameter can be obtained 

2.71 < a < 2.8 . (76) 

Therefore, when dilaton freezes after inflation, the best range for the coupling /($) 
can be determined from eq.(76). Considering / = /qO^""^, a = c t^^"^ (c is a constant) , 
/(^i?) = 1 and tji ~ ^ — , we have 

/o = c* {2H,^fY-^ . (77) 

Obtaining the value for /o and having the range of a, we can determine the dilaton 
coupling in this case. 



From Eq.(31) and the consistency condition that leads to v' — u — 1/2 = a — 1/2, 
we get 

^ = (-) 

and from Eq.(70), we have ^ ~ ^ that is close to the value ^, as derived by Turner and 
Widrow in [5]. Since a and ^ parameters are specified, the coupling to electromagnetic 
field for generating the lO'^^G magnetic field is determined. 

Conventionally, ^ is a constant parameter [5] and the term ^RA"^ will become very 
small after infiation since Ricci scalar in FRW metric is: R — °°+° +-^ and it will be 

a'' 

very small in the radiation and matter dominated epoch, hence it will be removed. 



6 Dilaton decay 

By now, we have assumed that the dilaton evolution ends after inflation. Let us now 
consider a more realistic situation (like in [13]) in which the dilaton continues its evolu- 
tion and reaches its minimum potential, then it begins to oscillate and then decays (the 
dilaton evolves from $«(< 0) to its minimum in $ = and here, f — 1). 
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6.1 Dilaton evolution after inflationary stage 

During coherent oscillation, the energy density p$, of the dilaton field evolves as 
which is slower than the energy density of radiation produced by the infiaton field via 
its decay, p^. If the condition p$ < holds until dilaton decays, the entropy per co- 
moving volume remains approximately unchanged, but if it changes to > Pr, a large 
amount of entropy will be produced [12]. If this happens, the magnetic energy density 
will be diluted by entropy production. 

We consider tosc — tT^'^ (w is dilaton mass) as the time when coherent oscillations 
begin. When t > tji, the universe is radiation dominated and thus, a(t) = aRit/tfjy^"^ . 
Before coherent oscillation, tn < t < tosc, the dilaton field evolves with its exponential 
potential. 

The amplitude of the dilaton field, ^r, is relatively small at the end of inflation. 
This is a consequence of the condition we imposed on the magnetic field energy density 
(consistency condition). When tR < t < tosc, the magnetic energy density increases by 
the dilaton evolution. The magnetic energy density should be smaller than the dilaton 
one on all scales so that it does not affect the dilaton evolution. As will be shown, larger 
values of Ak|$/j| result in larger values of the magnetic energy density and, thus, we get 
an upper limit for Ak|$r| from here. By solving the dilaton equation of motion with 
exponential potential (eq.(lO)) numerically, we get a range for Afi;|$j^| that satisfies these 
conditions [13]. The result is that \^r\ with A ~ should be of the order 1/k. 

With these conditions and the fact that V cx'p{—Xk^r)/ Pcj^ ~ ^ 1, the dilaton 
energy density will not be larger than the radiation energy density and we will, therefore, 
have entropy production only in the coherent oscillation epoch, t ~ tosc- By considering 
the slow roll condition, we have p$ ~ ^[^] ioi the dilaton field. The dilaton gets 
its minimum of potential at tosc and thus, p$ V. The minimum arises when other 
contributions from gaugino condensation enter into the to dilaton potential and generate 
a minimum [13]. If the universe stays radiation-dominated while the dilaton decays, 
rj\ then 

Pr^'^^H**) > (79) 

or 

P4^r > p4tosc)[^r (80) 

Or 0,[tosc) 



^>(^V^V. (81) 



Since, p$ V, one gets 

V \ m J \ m 

where, we have used r$ ~ m(m/MpiY and tR ~ \Hinf. Entropy per comoving volume 
remains constant here [12]. Thus, the necessary condition for entropy production is 



V \ m J \ m ' 
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Suppose the dilaton energy density equals radiation energy density ai t — tc 

pA^]-' - P^itosc)[^r . (83) 

We, therefore, have 

/n.x2/ 1 \4 

. (84) 



V \VJ \2H, 



''it ''OSC 



'in/. 



where pq,{tosc) ~ and t^^c — ^ arc assumed. After tc, dominates over pn and 

_1 _1 2 

universe becomes matter-dominated and thus a{t) = ac{t/tc)'^'^ = aut^^tc "ts that sec- 
ond equahty comes from continuity condition of a{t) aX t — Ir. 

The entropy per comoving volume is S = a^{p + p) /T where p, p and T are energy 
density, pressure, and temperature, respectively, in equilibrium [12] and we have 

S'/' = S^'/' + ^P^gad^^l^/^r^ A^] exp[-r^(r - Q]dT , 

O 45 Jtc CLc 

poo ^ 

s,'/'[i + r^t-'/' / {u + Q'/'eM-^^u)du] « s^'/'ii + rQitj'^)-'/'] . (85) 

Jo 6 

where Sc is the entropy per comoving volume at tc and < 51* > is the average of 
(number of degrees of freedom ) for decay duration. In the second approximation, 
we have u = t — tc and we consider the limit u — > 00. In addition, we have used 
Pritc) — P^{tc) and the following equation 

S^[lC-^y/'a%f^. (86) 

that in general, gives the relation between the radiation energy density and the entropy 
per comoving volume. Prom r$ ~ m{m/Mp{f and Eqs.(84) and (85), we find that the 
ratio of entropy per comoving volume after decay to that of before decay is written as: 




By entropy production, the universe expands more rapidly to cancel this entropy 
production. It follows from this observation and also from the relation p,. oc a^^S"^/^, 
that ^ 

(A^)4/3 =H\ (88) 
V / 

where is the present scale factor when we have entropy production. 
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Finally, we investigate the effect of dilaton decay on the energy density of present 
large scale magnetic fields. Again, we assume that the universe immediately becomes 
highly conductive after reheating. From Eqs.(52),(87), and (88), we find the ratio of 
magnetic field energy density at the present time, when dilaton continues evolving after 
reheating, p', to the situation that dilaton freezes in the reheating, p, is 

^ ^ r\tn)(ASr'/' , (89) 
P 

We define 

M (90) 

Pct> 

Since we assumed ptj, >> p$, then uj « 1 which we consider QjS (jJ 
From Eq.(89) and ^ ~ cuexp [Xk^r), we arrive at 

^^f-\tn)eMl XK^n][w{^n^)r/' . (91) 

As can be seen from RHS of the above equation, it is the dilaton coupling that 
increases p' relative to p. 

6.2 Magnetic fields with dilaton decay 

Considering / = /oa^"^^, a = c t^^"^ (c is a constant), /(tosc) = 1 and tosc ~ we have 

/o = c* m"-^ . (92) 

since tR=jj^, 

fitR)^m^-' (2//,„/)i-V (93) 

We estimate from Eqs.(52), (91) and (93), the present strength of magnetic fields for 
the case in which dilaton continues evolving after reheating which leads to the following 
relation: 

Log{B'{to)) = ^[{a - l)Log{H,af) + c(2 a - 6) + 2 Log{T{a - 1/2) 
+0.6( a - 1/2)) - 23.96 + Xk\^r\ + (a - l){Log{2Hinf) - Log{m)) + A5] . (94) 

c is equal to 16.041 and 17.041 on 1 Mpc and 10 Mpc scale, respectively. It can 
be seen that a stronger magnetic field results from a larger value of A/t|$i?|. As said 
before, Ak|$^| is of the order 1 and any bigger value violates the consistency condition. 
Therefore, we use the maximum possible value for Xh\(^r\{= 1) and plot the graph for 
Log{B') versus a (fig(3)). The graphs of magnetic field strengths are on 10 Mpc scale 
in "A" and "B" and on 1 Mpc scale in "C" and "D". We have entropy production in all 
the cases. The entropy production in "B", "D", and "F" is 4.4 x 10^ while it is about 
2 X 10^ in "A", "C", and "E". 



18 



Comparing these graphs with those of Fig(2), one can find that in this case smaller 
values of a give the desired magnetic fields (10~^^, 10"^'', 10"^° and 10~^G) and we don't 
have consistency condition violation in any of the field strengths. 

The graphs "E" and "F" represent a strength of lO"-*^^ (which is the best choice ) 
with Hinj — 10^^ Gev and Hinf — 10^^ Gev^ respectively. 

From graphs "E" and "F", the best range for a to give the acceptable strength is 
2.215 < a < 2.441 and when the value of ^ is obtained from Eq.(70), the couplings are 
completely determined. 

7 conclusion 

In this work, a mechanism is introduced for amplifying quantum fluctuations in the in- 
flationary universe and for producing large scale magnetic fields. This is accomplished 
by breaking the conformal invariancc of electromagnetic theory in the early universe by 
coupling the gravity (curvature of spacetime) and a scalar field (dilaton) to it. These 
couphngs have been studied separately, before. Here, the more realistic problem to the 
effect that both of these couphngs exist together is considered. 

Two different situations have been discussed for the evolution of the dilaton scalar 
field. In the first situation, whereby the dilaton is assumed to freeze at the end of infia- 
tion, the parameters of the model are determined for the seed magnetic field that gives 
the presently observed strengths by the effect of amplification mechanisms. By consid- 
ering the strength of 10~^^ gauss as the best value for the seed magnetic field to give 
presently observed field, we fixed our coupling parameters as : 2.71 < a < 2.8 (we use 
"2.71" that corresponds to the upper limit of i^m/ which is more realistic) and ^ ~ 

In the second situation, which assumes that the dilaton continues its evolution after 
the infiation, a large amount of entropy can be produced that dilutes the energy density 
of the magnetic fields produced. Here, we expect less values for a since our coupling 
nears unity in a larger time than the previous case. The range for variations of the pa- 
rameters that give the desired magnetic field ( 10~^^ gauss ) includes 2.215 < a. < 2.441 
and ^ ^ which is closer to the value derived in [5] . 
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Figure 1: Log{Hinf) relative to v', for magnetic fields with strengths 10 10 10 
on 1 MPc scale and 10~^G on 10 MPc scale, at present. 
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Figure 2: Log{Hinf) relative to a, for magnetic fields with strengths 10 10 10 
on 1 MPc scale and 10~^G on 10 MPc scale, at present 
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Figure 3: Log{B'{To)) relative to a, for Hinf = 2.4 x lO^^GeV, m = 2.4 x lO^^GeV in 
"A", "C" and "E", for Hi^f = lO^^GeV, m = lO^^GeV in "B", "D" and "F". Graphs 
"A" and "B" are on 10 Mpc and others are on 1 Mpc. Graphs "E" and "F" are depicted 
to get the accurate value of a. 
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